Large variations in the hole spin splitting of quantum- wire subband edges 
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We study Zeeman splitting of zone-center subband edges in a cylindrical hole wire subject to a magnetic field 
parallel to its axis. The (?-factor turns out to fluctuate strongly as a function of wire-subband index, assuming 
values that differ substantially from those found in higher-dimensional systems. We analyze the spin properties 
of hole-wire states using invariants of the spin-3/2 density matrix and find a strong correlation between g- 
factor value and the profile of hole-spin polarization density. Our results suggest possibilities for confinement 
engineering of hole spin splittings. 
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Spin splitting of charge carriers in semiconductors has been 
a focus of recent research interest, partly because it may form 
the basis for the new paradigms of spin-based electronics and 
quantum information processing. 1 Besides such possible ap- 
plications, intriguing fundamental-science questions motivate 
the study of charge carriers' spin properties. In particular, 
the quantum-mechanical coupling between spin and orbital 
degrees of freedom enables a host of, sometimes counter- 
intuitive, mechanisms for manipulating spins in nanostruc- 
tures. 2 As states in the valence band of a typical semicon- 
ductor are subject to a strong spin-orbit coupling, hole spin 
splittings will be highly tunable by geometrical and quantum- 
confinement effects. Large anisotropics of hole g-factors in 
quantum wells, 3 point contacts, 4 quantum dots, 5 and local- 
ized acceptor states 6 provide pertinent examples. The origin 
of such peculiar hole-spin properties can be traced to the fact 
that quasiparticles from the top-most valence band are charac- 
terized by total angular momentum (spin) 3/2. (Conduction- 
band electrons are spin- 1/2 particles like electrons in vacuum.) 
Although spin is an intrinsically quantum degree of freedom, 
it has been possible to rationalize spin- 1/2 physics largely in 
terms of a magnetic-dipole analogy. This classical-physics- 
inspired vehicle for our understanding succeeds because any 
spin- 1/2 density matrix is fully characterized 7 by the (trivial) 
particle density and a dipole moment associated with spin- 
polarization. In contrast, the spin-3/2 density matrix has two 
more invariants, a quadrupole and an octupole moment. 8 As a 
surprising implication, magnetic fields do not always induce a 
spin polarization in two-dimensional (2D) hole systems. 9 

Our theoretical study presented here reveals the drastic in- 
fluence of a quantum-wire confinement on spin-3/2 physics. 
Figure Q] shows universal results for the Lande <?-factors of 
low-lying subband edges in cylindrical quantum wires subject 
to a parallel magnetic field. Two surprising features are appar- 
ent. First, g* is seen to assume unusual values. Considering 
the wire axis to be a natural quantization axis for hole spin and 
remembering that spin-3/2 projection eigenvalues are ±3/2 
and ±1/2, we would expect to find only 6k and 2k as possible 
g-factors. (Here k is the hole g-factor in the bulk material.) 10 
In Fig.Q] quite different numbers are found. Second, a strong 
variation of g-factor values between different subbands is ob- 
served. These features contrast also with situations found in 
2D hole systems where confinement results in the suppression 



of field-induced spin polarization 3 and g-factors close to and 

4k emerge d 1 ' 12 ' 13 

In the following, we show how the anomalous behavior 
of wire-subband spin splittings is rooted in peculiar patterns 
of hole-spin dipole density that arise from the interplay be- 
tween quasi-onedimensional (ID) confinement and spin-orbit 
coupling in the valence band. Our results are experimentally 
testable, as Zeeman splitting of quasi- ID subband edges can 
be measured directly in transport experiments. 14 P-type ver- 
sions of recently fabricated highly symmetric semiconductor 
nanowires 15 should be particularly well-suited for observing 
effects discussed in this work. Our study is also relevant for 
understanding spin properties of quasi- ID hole systems real- 
ized in 2D semiconductor heterostructure a 4 ' 16 ' 17 and provides 
physical insight for the interpretation of recent numerical re- 
sults^jo 
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FIG. 1: (Color online) Lande factors g* of low-lying subband edges 
in a cylindrical hole quantum wire (in units of the bulk hole <?-factor 
k). The abscissa shows subband energies measured from the valence- 
band top in units of Er = — 7ift 2 /(2mo-R 2 ). 71 is a materials- 
dependent (Luttinger) parameter, mo the electron mass in vacuum, 
and R the wire radius. The calculation is based on the spherical 
Luttinger model 10 for the top valence band and assumed the GaAs 
value (0.37) for the relative spin-orbit coupling strength 75/71- 
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The remainder of this article is organized as follows. We 
first describe the theoretical formalism for calculating hole- 
wire subband energies and g-factors. After that, the character- 
ization of spin-3/2 density matrices in terms of multipole in- 
variants is introduced and results for the cylindrical-wire sub- 
band edges are presented. The emerging correlation between 
hole spin-dipole density and g-factor will be discussed. In 
closing, we comment on the effect of deviations from certain 
idealities assumed in our theoretical description. 

Cylindrical hole quantum wire in a parallel magnetic field. 
We adopt the Luttinger model 10 to describe states in the top- 
most valence band of a bulk semiconductor. The correspond- 
ing Hamiltonian operates in the 4x4 space spanned by eigen- 
states of spin-3/2 projection on an arbitrary axis with quantum 
numbers ±3/2 (heavy holes, HH) and ±1/2 (light holes, LH). 
It reads in the spherical approximationiSi 2 ! 
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We denote linear orbital hole momentum by p, J is the vec- 
tor of spin-3/2 matrixes, and too the vacuum electron mass. 
Energies are measured from the valence-band edge. The ratio 
7 = 7 s /7i of Luttinger band-structure parameters quantifies 
the relative strength of spin-orbit coupling in the valence band. 

We consider a wire that is parallel to the z axis and sub- 
ject to Zeeman splitting due to a magnetic field pointing in 
the same direction. Thus it is useful to choose the represen- 
tation where J z is diagonal. As we are focused on finding 
hole-wire subband edges at the zone center, 22 we set p z = 0. 
Orbital magnetic-field effects are neglected. 23 To determine 
bulk hole states that can be superimposed to find cylindrical- 
wire subband states, we use polar coordinates r, <p and the 
wave-function ansatS^ 



that can be superimposed to find wire-subband bound states: 

i>kma+{r) = (a Q+ J m (fcr),0, J m+2 (fcr),0) T , (4a) 
i>kma-(r) = (0, J_ m _ 2 (fcr),0, a Q _J_ m (fcr)) T . (4b) 

We omit the polar-angle (<ys)-dependentpart since the wire po- 
tential is cylindrically symmetric. In Eq. (0), a = ±1 labels 
eigenstates within subspaces spanned by HH and LH basis 
states having spin projection ±3/2 and =Fl/2 (for ^kma±, re- 
spectively). The coefficients a aa are given by 
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l-2aJl + C ka +( 2 ka + 2Cka 



(5) 



where we used the abbreviation (ka — oneB /(-f s hk 2 ), and 
— e is the electron charge. The corresponding eigenenergies 
are, like the a aa , independent of orbital angular momentum 
quantum number m: 
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These results generalize previous analytical treatments 2 ^ of 
cylindrical hole wires to the case with finite Zeeman splitting. 

Inspection of the bulk dispersions in a finite magnetic field 
(Fig. |2]i reveals an interesting cross-over behavior. At small 
k, i.e., small kinetic energy for motion perpendicular to the 
magnetic-field direction, a Zeeman splitting of HH and LH 
states with Lande factor in accordance with naive expectation 
(Shh = ® K and <j£ H = 2k) is found. Evidently, these states 
have their spin closely aligned with the z axis. In contrast, 
at large fc, the less-dispersive branches show no Zeeman split- 
ting, while that of the more dispersive ones is characterized by 
g-factor An. The latter situation is reminiscent of the Zeeman 
effect in 2D hole systems for an in-plane magnetic field. 3 
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J m {x) is a Bessel function, and a m , . . . , d m are constants. 
Uniqueness of the hole wave function requires m to be an in- 
teger. Inserting Eq. (f2]i into the stationary Schrodinger equa- 
tion simplifies the Hamiltonian for hole motion transverse to 
the z axis, which now reads 24 
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(3) 

Here J± = (J x ± iJ y )/ \[2 are ladder operators for the z-axis 
projection of spin. 

Zeeman splitting due to a magnetic field with magnitude B 
applied parallel to the z axis is described by Tiz = 2k[IbBJ z , 
which needs to be added to Eqs. (Q]) and J3). \lb is the Bohr 
magneton, k the bulk hole g-factor, and we neglected the small 
anisotropic part of Zeeman splitting in the bulk valence band. 

Diagonalization of 7i\^ + Hz yields bulk-hole eigenstates 




FIG. 2: (Color online) Energy dispersion E aa (k) of bulk-hole states 
that can be superimposed to form wire-subband bound states. Here 
R is an arbitrary length scale that will later be identified with the wire 
radius. The energy unit is En = —71ft 2 / (2mo-R 2 ), and a magnetic 
field B = 10 • 7ift/(2Kei? 2 ) has been applied parallel to the z di- 
rection. The GaAs value 7 = 0.37 was assumed for the calculation. 
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TABLE I: Energies and effective (/-factors for the lowest ten quasi- 
1D subband edges in a cylindrical hole wire. Units and parameters 
are the same as given in the caption of Fig.Q] 
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Hole-wire subband bound states are found by superimpos- 
ing spinors ipk atT maa(r) with fixed m, a and E aa (k aa ) = E 
to satisfy a hard-wall boundary condition at r = R. The re- 
sulting secular equation reads 

a + J m (k + R)J m+2 (k^R) ~ a-J m (k^R)J m+2 {k + R) = 0, 

(7) 

with conserved index a suppressed for the sake of brevity. 
Solving Eq. (0 for E yields the subband energies E nmrT (B) 
for a cylindrical hole wire subject to a parallel magnetic field 
of magnitude B. In general, these are mixtures of the bulk 
states "4>km±rr- Exceptions are states with m = — 1, for which 
Eq. (Q specializes to the condition Ji(k±R) = that can 
be satisfied by the individual states tpkm+a and ijjkm-a- At 
finite B, levels E nm ± are Zeeman-split. We calculate the cor- 
responding g-factors using2£ 



= lim 



E n m+ (B) — E nm _ (B) 



p B B 



(8) 



Table U summarizes results obtained for the ten lowest wire 
subbands when the GaAs value for 7 is assumed. The strong 
variation of g* for different wire levels is illustrated in Fig.Q] 

The first subband has g* close to 2k, as expected for a LH 
state with spin polarization parallel to the wire axis— The 
second level belongs to the above-mentioned special class of 
states with m = — 1, which are pure bulk states and therefore 
have g-factor equal to either or 4k. To obtain a fuller under- 
standing of how various effective g-factor values for hole-wire 
subbands edges emerge, we analyze the corresponding bound 
states in terms of invariants for the spin-3/2 density matrix. 

Spin polarization of hole-wire bound states. The peculiar 
physical properties of hole quantum wires can be attributed, 
in a very general way, to finite HH-LH mixing present even 
at the subband edges. 27 A quantitative characterization of the 
latter is typically attempted by considering expectation values 
for spin-3/2 projections on a fixed axis i 18,20 ' 27 In particular for 
a wire geometry, the usefulness of such an approach can be 
limited by ambiguities in the choice of a suitable projection 
axis. Hence, a characterization of hole states in terms of scalar 
invariants would be much more meaningful. In the following, 
we provide such an analysis in terms of a multipole expansion 
of spin-3/2 density matrices. 8 

As hole states are four-spinors, there exist four scalars 
po, . . . , ps to characterize their spin-3/2 density matrix. Here 
we adopted the notation of Ref.8, where po is proportional to 
the hole (charge) density, p\ is a dipole moment related to the 
hole-spin polarization, the quadrupole moment p 2 quantifies 
HH-LH mixing, and p% is an octupole moment. It is an in- 
triguing feature of spin-3/2 physics that magnetic fields can 



induce a substantial octupole moment instead of a spin po- 
larization (dipole moment) in 2D hole systems. 8 Here we ob- 
serve an analogous property of hole spin for individual states 
at the quasi- ID subband edges and relate our findings to the 
measurable g-factors of these states. 

Within our model for a hole quantum wire, states at the sub- 
band edge are superpositions of HH and LH amplitudes cor- 
responding to J z projection cr3/2 and — a 1/2, respectively. 
[See Eqs. ©.] For these wave functions, relations exist be- 
tween the multipole invariants of the spin-3/2 density matrix: 



pi = 2p\ =pl + p 2 3 



(9) 



The left equality in Eq. (0 quantifies the HH-LH mixing that 
is present in subband-edge states. The right equality in the 
same equation implies that, for each of these states, the mag- 
nitudes of spin polarization and octupole moment are com- 
plementary. In particular, a state with no spin polarization 
will have a maximum octupole moment and vice versa. We 
are thus able to fully characterize the spin properties of each 
subband-edge bound state by considering only p\j p\. When 
derived from the radial part of the wire-subband wave func- 
tion, this quantity provides a measure of local hole spin polar- 
ization. Figure [3] shows results for the five lowest cylindrical- 
wire subband edges. The (mainly) LH character of the lowest 
subband is apparent (p\j ' p\ is almost constant at 0.2), as is the 
vanishing of spin polarization for the next subband. Higher 
subbands show an evolving mixture of HH/LH character in 
their spin-polarization profile. 

A comparison of spin-polarization profiles (Fig. [3} and ef- 
fective g-factors (Table [Ji for each level yields a consistent 
picture. States with mainly LH (HH) character, i.e., those hav- 
ing p\j p\ ~ 1/5 (9/5) can be associated with a definite spin 
polarization parallel to the wire axis. Their corresponding g* 
is close to the expected value 2k {6k). The lowest subband- 
edge state is an example for such a polarized (LH) state. Other 
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FIG. 3: (Color online) Radial profile of the (normalized squared 
magnitude of) local hole-spin polarisation, calculated for the five 
lowest wire-subband-edge bound states. Pure HH (LH) states, de- 
fined by spin projection parallel to the z axis (i.e., the wire axis), 
have pllpl = 9/5(1/5). 
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states exhibit a mixed spin-polarization profile, i.e., have both 
HH and LH polarizations present in different regions across 
the wire. The associated (7-factor values can vary widely. Fi- 
nally, states exist with vanishing hole-spin polarization over 
(most of) the wire cross-section. By virtue of Eq. (O, these 
states have a large spin-3/2 octupole moment. The second- 
lowest subband is an example for this class of states, whose 
(7-factor turns out to be (close to) zero. The remarkable cor- 
respondence between local spin-polarization profiles and g- 
factor values for hole-wire levels is another example for in- 
triguing spin physics emerging in nanostructures with strong 
spin-orbit coupling. The possibility to access individual prop- 
erties of quasi- ID subband edges in transport experiments^ 4 - 
should enable the detailed study of hole states with large and 
small spin polarizations and spin octupole moments, respec- 
tively. 

Principal caveats. To gauge the relevance of our results 
for real experiments, we need to discuss a few idealizations 
and approximations that are implicit in our model. Most im- 
portantly, we neglected orbital magnetic-field effects, band 
warping, and coupling to the conduction and split-off valence 



bands. In addition, realistic wires will show deviations from 
perfect cylindrical symmetry and usually have a softer than 
hard-wall confinement. 

In principle, any physical property affecting HH-LH mix- 
ing will quantitatively change g*. It can be expected, how- 
ever, that effects of band warping and remote bands are small, 
as they turned out to be for 2D hole systems. 8 Cross-sectional 
shape of a hole wire generally affects its spin splitting. 26 How- 
ever, a recent numerical study 20 found results for square CdTe 
wires that agree closely with those presented in this work. A 
detailed comparison with experiment will need to address or- 
bital effects due to the typically not-so-small magnetic fields 
used to measure ff* i 4 ' 14 i 23 Interpretation of data obtained in 
hole point contacts 4 also requires a better understanding of 
the transition region between 2D HH contacts and a constric- 
tion, and a proper treatment of the strong quantum-well and 
soft lateral confinements. 
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